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Abstract 

Exoplanets with lower equilibrium temperatures than Earth and primordial hy- 
drogen atmospheres that evaporate after formation should pass through transient 
periods where oceans can form on their surfaces, as liquid water can form below 
a few thousand bar pressure and H2-H2 collision-induced absorption provides sig- 
nificant greenhouse warming. The duration of the transient period depends on the 
planet size, starting H2 inventory and star type, with the longest periods typically 
occurring for planets around M-class stars. As pre-biotic compounds readily form in 
the reducing chemistry of hydrogen-rich atmospheres, conditions on these planets 
could be favourable to the emergence of life. The ultimate fate of any emergent or- 
ganisms under such conditions would depend on their ability to adapt to (or modify) 
their gradually cooling environment. 
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1 Introduction 



Recent developments in exoplanet characterization [e.g., Tinetti et al. , 2007 



Knutson et al. 2007, Bean et al. 2010 have raised the possibility that in 



the relatively near future, we may be able to search for biosignatures in the 
atmospheres of terrestrial-mass planets. However, as Earth is the only current 
example of a life-supporting world that we know of, there is still no consensus 
as to the form such biosignatures will take. More fundamentally, there are 
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still major unknowns as to the prevalence of conditions necessary to create 
and sustain life on other planets in general. 



To date, much of the discussion has centered on the concept of the habitabil- 
ity zone - the orbital range around a star in which temperatures are in the 
right range for a planet to have surface liquid water. The standard definition 
of this region Kasting et al. , 1993 is based on the assumption that a habit- 
able planet will have a climate dominated by the greenhouse warming of CO2 
and H 2 0, with temperature regulation provided by the carbonate-silicate cy- 



cle via plate tectonics. Alternative hypotheses exist: for example, Stevenson 



1999 proposed that if they retain primordial H2-He envelopes, planetary em- 



bryos ejected from their systems during planet formation could be kept warm 
enough for surface liquid water even as they drift through interstellar space. 
Following this, Pierrehumbert and Gaidos 2011 noted that orbitally bound 
exoplanets far from their host stars might also be habitable to life if they 
form with relatively thin hydrogen envelopes and manage to preserve them 
against escape. Both these studies suggest that a definition of the habitable 
zone based on H2O-CO2 warming only may be too narrow. However, the per- 
centage of planets left with exactly the right atmospheric pressures to allow 
habitable temperatures after the early stages of atmospheric erosion is likely 
to be small, and the observation of rocky planets that are free-floating or in 
distant orbits around other stars is technically challenging. 



Here the impact of hydrogen greenhouse warming is considered for exoplanets 
that are closer-in and hence may lose their hydrogen atmospheres entirely. 
Major uncertainties remain in planetary formation theory, but it is generally 
understood that planets of between 1 and 20 Earth masses can form with 
primordial hydrogen envelopes either due to capture of the gas from the so- 
lar nebula Rafikov, 2006 or outgassing during accretion Elkins-Tanton and 



Seager, 2008 . In the case of Earth and Venus, these envelopes may have been 



lost to space extremely rapidly, while as noted by Stevenson 1999 and Pier 



rehumbert and Gaidos 2011 , planets that are in distant orbits or are ejected 
to interstellar space during formation can retain them permanently. In the 
intermediate situations where hydrogen atmospheres persist some time after 
formation but are eventually removed, the planet's climate will clearly evolve 
over time. Here I argue that during this evolution, transient periods will occur 
where water exists as a liquid on the planet's surface, and organic molecules 
may form in the atmosphere via photochemistry. 



In the following section, the greenhouse effect of hydrogen atmospheres is 
discussed and an empirical expression for the pressure necessary to maintain 
surface temperatures in the correct range for liquid water is derived. Next, a 
simple assessment of atmospheric loss due to thermal escape is presented, and 
transient habitable times as a function of planet mass, orbital distance and 
starting hydrogen inventory are estimated. Finally, the effects of processes not 
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included in the model are considered, along with the potential implications of 
the results. 



2 Analysis 



2. 1 Climate 



In atmospheres dominated by hydrogen, the main opacity source in the in- 
frared is H2-H2 collision-induced absorption (CIA). Other gases can also con- 
tribute to infrared opacity and hence greenhouse warming, but only if they 
are present in sufficient amounts in the mid- to high atmosphere to absorb in 
the spectral regions where hydrogen becomes transparent. Helium can warm 
through He-H 2 or He-He CIA, but its effects should be secondary to those of 
hydrogen. In addition, it will not be present in large quantities if the atmo- 



sphere is created by accretion outgassing Elkins-Tanton and Seager, 2008 



For temperatures of a few hundred Kelvin or less, the dominant equilibrium 
O-C-N species will be H 2 0, CH 4 and NH 3 if hydrogen is the majority gas, as 
in the atmospheres of the Solar System gas giants [e.g. Moses et al. 201 1 



The primary effect of these gases should be to increase greenhouse warming, 
although they can also influence the climate via cloud and haze formation (see 
Discussion). For simplicity, the climatic effect of gases other than H 2 is ne- 
glected in the main part of this analysis, although some calculations involving 
H 2 -H 2 are also discussed. 



To allow rapid assessment of climate for a range of cases, one-dimensional 
radiative-convective calculations were used to produce an analytical expres- 
sion for the globally averaged outgoing longwave radiation (OLR) F out and 
shortwave albedo A assuming pure H 2 atmospheres. The model used was sim- 

2010a|b| : there 



ilar to that described in Wordsworth et al. 



bands in the longwave and 36 in the shortwave, and AD Leo and solar data 
|Segura et al. , 2005 were used for M-class and G-class stellar spectra, respec- 
tively. A two-stream scheme Toon et al. , 1989 was used to compute fluxes, 



and the radiative effects of Rayleigh scattering in the visible were taken into 



account using the H 2 refractive index from Hecht 2001 



Simulations were performed with specific heat capacity c p = 14.31 kJ kg 1 
K -1 Lide, 2000, value at 300 K] and molar mass fin 2 = 2.02 g to calculate 



F out as a function of surface pressure p s and gravity g. Hydrogen CIA opaci- 



ties were calculated from tabulated spectral absorption data Borysow , 2002 



The vertical temperature profile was assumed to follow the dry ideal adiabat 
T(z) = T s — (g/c p )z, with the transition to an isothermal region ('strato- 
sphere') at low pressures assumed to occur at 100 K (Figure [T]). Comparison 
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with a time-stepping version of the model and simulations where the isother- 
mal transition was neglected entirely showed that the stratospheric temper- 
ature did not have a significant influence on the value of F out , in agreement 
with previous studies |Pierrehumbert and Gaidos , 2011 . 



We are most interested in greenhouse warming for surface temperatures in the 
range T min <T S < T max , where {T min ,T max } = {273,333} K are conservative 
temperature limits for habitability to Earth-like microbial life. The variation 
of the effective OLR F out /aT^ with temperature was found to be slow, so cal- 
culations were performed at T s = 300 K only and the temperature dependence 
of the effective OLR was neglected. Given this simplification, the expression 

F out /aTt = ap- b g b J 2 (1) 



was calculated from the numerical results within the limits 0.5 < p s < 10 bar 
and 10 < g < 30 ms -2 using a logarithmic least squares fit (see Figure [2|. 
Here p s is in bars, a = 0.279, b = 1.123, g* = g/g re f and g re f = 20.0 ms" 2 . 

To determine the potential error in the analysis due to the assumption of pure 
H2 atmospheres, a few additional calculations were performed including the 
presence of water vapour. In this case, saturation of H2O vapour in the atmo- 
sphere was assumed, and the lapse rate was calculated with H 2 assumed to 
be an ideal gas. The effect of H 2 on the radiative transfer was included in 



a similar way to in previous studies [e.g. Wordsworth et al. , 2010b . High res- 



olution spectra were computed from the HITRAN database iRothman et al 



2009 on a 6 x 8 x 8 temperature, pressure and H2O volume mixing ratio grid 
with values T = {100, 150, . . . , 350} K, p = {10~ 2 , 10"\ . . . , 10 5 } mbar and 
Qh 2 o — {10 -7 , 10~ 6 , . . . , 10°}, respectively, before conversion to a correlated-A; 
format for use in the radiative-convective model. 

Figure [T] shows the temperature profiles vs. pressure for pure H 2 and H2-H2O 
atmospheres with constant stratospheric temperatures of 100 K. Figure [2] 
shows the resulting effective OLR calculated using the radiative-convective 
model for a planet with g = g re f. As can be seen, water vapour changes the 
lapse rate in the low atmosphere, primarily via latent heat release, but its 
effects decrease rapidly as the total pressure increases beyond a few fractions 
of a bar (the vapour pressure of H2O is ~ 0.035 bar at 300 K). These lapse 
rate changes were more than compensated for climatically by increased in- 
frared absorption - as shown by Figure |2j water vapour generally decreased 
the effective OLR and hence increased the greenhouse effect. Nonetheless, its 
effects were not significant for pu2 values above about 0.4 bar. Hence equation 
(II]) was judged an acceptable approximation for use in the following analysis. 

In the visible, the radiative transfer in pure H 2 atmospheres is dominated by 
a combination of CIA and Rayleigh scattering. For surface pressures greater 
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than around 20 bar, the fraction of starlight reaching the surface is small, and 
calculation of the surface temperature becomes difficult without knowledge of 
the planet's internal heat flux. In this analysis, all climate calculations were 
restricted to atmospheric pressures of 10 bar or less, for which it was sufficient 
to determine the planet's bond albedo in the visible. A surface albedo of 0.2 
(typical for a rocky planet) was assumed. Albedos A ~ 0.25+0.151ogio(p s [bar]) 
and ~ 0.15 were found for G-class and M-class stars, respectively. The latter 
values were lower and approximately constant because the effects of Rayleigh 
scattering depend strongly on wavelength, and M-class stars have redder spec- 
tra. The variation of albedo with surface gravity was found to be small (under 
0.025 at T s = 300 K) and hence ignored. The effect of water vapour on albedo 
due to absorption was slightly greater (maximum decrease of ~ 0.1), but given 
the fact that it tended to increase warming and the inherent uncertainties due 
to haze and cloud cover, it was also neglected for simplicity. Finally, the re- 
duction in visible flux early in the lifetime of G-class stars Gough, 1981 was 
taken into account, with the G-star stellar luminosity taken to be L = 0.7 X Q . 
In the M-class case, the value L = O.O12O5L (representative of the red dwarf 
GJ581; von Braun et al. 2011] ) was used. 



Planets with slowly escaping atmospheres will be in thermal equilibrium, 
Font = F in , if the loss timescale ti oss is much longer than the radiative adjust- 
ment timescale t ra d = c p p s T s / gaT* |Goody and Yung 1989 where o is the 



Stefan-Boltzmann constant, g is gravity and T e is the characteristic emission 
temperature. Hydrogen has a mean optical depth of 1 in the infrared at around 
0.2 bar pressure, so for e.g., a 20 bar atmosphere with surface temperature 
300 K, the emission temperature should be T e = T s (p e /p s ) R/ ^ CpflH ^ ~ 80 K, 
where R is the ideal gas constant. This yields t ra d = 12, 000 years for an 
Earth-gravity planet. The presence of a primordial steam atmosphere could 



complicate this picture [e.g. Valencia et al. , 2008 , but only for planets that 



lose their entire hydrogen inventory in the first few million years after forma- 
tion. In the following analysis this possibility is neglected and F out = F in is 
assumed. 

Writing Fi n = |(1 — A)(L/47id 2 ), where L is the stellar luminosity and d 
the orbital distance, and making the standard ID modelling assumption of 
efficient horizontal heat transport by the atmosphere, we find 



G \ ,1 /no / (ICfTf 



p. 



M^-n - Z ' * T, \ , (2) 



\v 2 g re f) \(l ~ A)L/lQ7rd 2 

(xM 0.23 T 3.56 rf 1.78 (1 _ A) -0.89 



2 Here we use T s rather than the smaller T e as an upper limit estimate of the mean 
atmospheric temperature. 
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Sotin et al. 



where the mass-radius relation r = (ar e /M§)M P = vM? from 
2007| is used to write g in terms of the planetary mass M. For simplicity 
the rocky planet values a = 1.0 and (3 = 0.274 are assumed throughout 
the calculation. The approximately constant albedo in the M-star case allows 
direct calculation of p s [M, d, T s ], while in the G-star case, iteration is necessary 
for a self-consistent calculation. For an example 5M e rocky planet at 2.5 AU 
around a Sun-like star with L = 0.7 X Q , iteration of (|2| predicts an albedo 
of 0.35 and 4.3 bars of H2 necessary to maintain a surface temperature of 
300 K. The saturation temperature of H2 at this pressure is less than 25 K 
so atmospheric condensation does not occur at any altitude. As 

have already noted, 



Lide 2000 



Stevenson 1999 and Pierrehumbert and Gaidos 2011 



hydrogen is an extremely potent, effectively incondensible greenhouse gas. Its 
main disadvantage in the context of achieving habitable temperatures is of 
course that usually, either far too much or too little of it is present. 



2.2 Atmospheric escape 



Unless hydrogen is continuously supplied via vulcanism or other surface pro- 
cesses at the same rate at which it escapes, habitable conditions due to H 2 -H 2 
warming will clearly be transient on exoplanets with escaping atmospheres. 
The habitable duration will depend primarily on the host star's extreme ul- 
traviolet (XUV) emission and the planet's size and orbital distance. Erosion 
due to the stellar wind and meteorite impacts can also play a role in hydrogen 
loss from rocky planets, but thermal / hydrodynamic escape processes should 
dominate in the early phase of a system's evolution when stellar XUV levels 



are high Tian et al. , 2005 



For an atmosphere undergoing thermal escape, the flux of hydrogen to space 
depends primarily on the thermal escape parameter A = r/H exo , where H exo = 
kBT eX o/rriHg is the scale height, ks is the Boltzmann constant, T exo is the 
temperature just below the exobase (exospheric temperature) and m H is the 
atomic mass of hydrogen^] The exospheric temperature is the most important 
unknown in A. It may depend on factors such as the XUV heating, coupling 
with the lower atmosphere via thermal conduction and wave propagation, and 
cooling to space due to radiatively active molecules / ions such as H3 [e.g., 



Koskinen et al. 


et al. 


2003 as 



2007 . Here T exo is calculated using the approach of Lammer 



T^xo — CFxuv/ 9 + T* 



(3) 



In the main part of this analysis efficient H2 dissociation is assumed, so that the 
dominant escaping species is atomic hydrogen. 
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where Fxuv is the stellar XUV flux, T min the temperature at the base of the 
thermosphere, s = 0.73 and C = 1.67 x 10 6 K s W 



m 3 s 



2 . In Lammer et al. 



2003), it was assumed following Bauer 1971 that T exo is primarily determined 



by a balance between XUV energy input and energy loss to the lower atmo- 
sphere due to thermal conduction, with the latter equal to K(T exo ) = K T^ XO 
ergs K _1 cm" 1 s _1 , where Kq = 16.4 for atomic hydrogen. This approach ne- 
glects the possible effects of radiative cooling and internal heat sources, which 
we discuss later. 

The temperature T min is included to account for conductive coupling between 
the thermosphere and the lower atmosphere. For all calculations assuming 
XUV-dominated escape, we set this parameter to the representative cold-trap 
temperature 100 K. However, the exact value of T m i n had a negligible influence 
on the results, as by the time the exospheric temperature drops to a few 
hundred Kelvin or lower, thermal escape rates are extremely small for planets 
of mass greater than ~ 1 mg. 

The stellar XUV flux depends primarily on the star type and age. Although 
XUV emissions may vary widely in a given stellar population (particularly 
for M-stars) the XUV flux tends to saturate at a value Fxuv ~ 10~ 31 L/47rd 2 
early in a star's lifetime, for a time txuv that ranges from around 100 million 
years (G-stars) to 1 billion years (M-stars) |Selsis et al. 2007 . After this it 

In 



Ribas et al., 2005 



decays relatively rapidly according to a power law 
most cases, the majority of atmospheric erosion therefore occurs early in a 
planet's lifetime. When other heating sources are negligible, an upper limit on 
atmospheric escape 



cFxuv 

4r 
eFxuv 



(4) 

[Pa s" 1 ] 



where e is an efficiency factor, can be derived. Equation Q assumes that the 
averaged XUV energy input to the atmosphere Fxuv I '4 (x the factor e) is used 
to overcome the gravitational potential energy of exospheric hydrogen atoms. 
For convenience when coupling with the climate equations, the hydrogen flux 
is converted to a loss rate of surface pressure per unit time. This makes use of 
the fact that surface pressure is equal to m# 2 <? times the area number density 
of H2 molecules, and the potential energy of one hydrogen atom is GMmn/r 
(with niH 2 ~ 2m h)- 

This simple energy balance approach gives a reasonable estimate of escape 
rates for values of A equal to or less than around 2.8, for which the escape can 



be regarded as hydrodynamic Volkov et al. , 2011 . For extremely small values 



of A, the outer layers of the atmosphere extend to several planetary radii and 
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the value of may be greater than that given in Q. Hydrodynamic escape in 
such atmospheres is a formidable problem involving three-dimensional effects 
and coupled physics/chemistry, which we do not attempt to tackle here. As will 
be seen, the frequency of transient habitable conditions is most likely limited 
by the total amount of hydrogen that can be removed from the atmosphere, 
so the use of (El) allows for a conservative estimate. 

For A > 2.8, the flux of hydrogen to space was estimated from the Jeans 
formula 

<fr = v p n exo f(X)m H g, [Pa s" 1 ] (5) 



where v p = (2/c B T e:ro /m#) 1//2 is the most probable particle speed, n exo = 1/ 
(aHH exo ) is the particle number density at the exobase, an = 7r(3.1 x 1CT 12 ) 2 m 2 
is the hydrogen atom collision cross- sectiorQand /(A) = (1 + A)^/ 2 ^ is 
the fraction of upward-moving particles with speed greater than the escape 
velocity. 

Figure [3] shows the total quantity of hydrogen that can be lost to space in 
time txuv as a function of planet mass and orbital distance for G- and M-class 
stars. Due to the relatively small radii of rocky planets and high exospheric 
temperatures under elevated XUV conditions, energy-limited loss governed by 
Q occurs instead of Jeans escape almost everywhere. As can be seen, more 
than 100 bars of hydrogen may be lost in G-class systems for orbits closer 
than 1.5 AU, while for M-class systems, the 100-bar boundary is 0.5-0.7 AU. 
Accounting for the slight difference in luminosities used, these results are in 



approximate agreement with those of Pierrehumbert and Gaidos 2011 , al- 



though they predict less escape than that shown in Figure [3] for more distant 



orbits. Pierrehumbert and Gaidos 2011 used a slightly more complex exo- 



sphere model that predicted the exospheric temperature numerically, but still 
neglected all significant radiative cooling effects. Onset of Jeans escape at 
lower values of d and M due to lower exospheric temperatures in their model 
is the most likely cause of the difference in the results. 

Given our simplifying assumption of constant saturated XUV flux throughout 
the time period txuv-, ©J @ and (|5]) can be combined to give an estimate of 
the duration of habitable conditions on the planet's surface 

At hab = [p(T max ) - p(T min ))/(f). (6) 



Figure [4] shows Athab in G- and M-class systems, assuming various initial 
H_2 inventories pq. In situations where there is too little atmospheric loss to 
reach p(T max ) after txuv, the results have not been plotted. As can be seen, 



Value at 1000 K, taken from 



Pierrehumbert 



2011 
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transient habitable conditions occur for many values of d, M and po- The 
habitable duration Athab varies widely, from a few thousand years for close-in 
planets to values approaching t X uv f° r more distant ones. In some cases, the 
amount of atmospheric erosion is just right to leave the planet with surface 
temperatures in the habitable range permanently (black regions in Fig. |1J). 
However, these cases represent only a small fraction of the total. 

The starting H2 pressure is a critical factor in the calculation. As Figure [4] 
shows, for large values of po only planets close to the star experience transient 



habitable periods. Pierrehumbert and Gaidos 2011 argue based on the results 



of Rafikov 2006 that initial H 2 pressures of 1-100 bar are possible, although 
it is likely that the exact amount in any situation will depend sensitively on 
the nature of accretion in the final stages of planet formation. As M-stars tend 
to have elevated XUV levels for longer periods than G-stars, planets around 
them can lose larger quantities of hydrogen, and hence have a greater chance 
of experiencing transient hydrogen-induced habitable periods given starting 
H2 pressures of greater than a few tens of bars. 

The influence of exospheric radiative cooling on escape were neglected in 
this analysis for simplicity, but it could nonetheless be important in some 
cases. The effect of cooling can be crudely investigated by fixing the ex- 
ospheric temperature at a constant value for XUV levels below the critical 
threshold Fxjjv ~ 0.18 W m~ 2 given by Koskinen et al.| |2007|. If we assume 



T exo = 3000 K for all planets outside a radius d cr u = y 10~ 3A L/4:irFxjjvi based 
on their Fig. 1, A < 2.8 then requires M < 1.16M©. For larger planetary 
masses, inefficient Jeans escape dominates and the amount of hydrogen re- 



moved from the atmosphere is significantly reduced. Unfortunately, Koskinen 



et al. 2007 examined only Jupiter-mass planets around Sun-like stars in their 
study, so we are unable to investigate the scaling of T exo with planetary mass 
and star type here. However, if Hjj~ cooling effects for terrestrial-mass planets 
are as strong as those for gas giants, transient habitable periods may be less 
common outside the critical radius than Fig. [4] suggests. 

For low-mass planets, thermal escape processes that do not depend on the 
stellar XUV flux may also contribute to H2 loss and hence eventual transient 
habitability. In the Solar System, the variation in mean exospheric tempera- 
tures for the gas giants is surprisingly small, with typical values in the range 
700 — 1000 K despite the large differences in solar XUV insolation. Internal 
processes such as heating due to magnetic field ion acceleration and gravity 



wave breaking are believed to dominate the energy budget Lunine, 1993 
If similar mechanisms were present for rocky exoplanets with hydrogen en- 
velopes, they would continue to lose gas to space even after the initial phase 
of saturated stellar XUV fluxes. Figure [5] shows the amount of hydrogen lost 
after 10 GYr due to atmospheric Jeans escape only, for planets of different 
masses and fixed T exo values. As can be seen, planets of Earth mass or lower 



9 



can continue to lose significant amounts of hydrogen even far from their host 
star if their exospheric temperatures remain above ~750 K. Such bodies will of 
course be difficult to observe, but the example serves to illustrate the variety 
of situations that may lead to transient habitable scenarios. 



3 Discussion 



In summary, if a planet with equilibrium temperature lower than that of Earth 
has multiple bars of hydrogen in its atmosphere early in its evolution (either 
due to direct capture from the nebula or outgassing after accretion) but subse- 
quently loses it due to escape, it will almost certainly pass through a transient 
state where liquid water can form on its surface. It has been shown that the 
duration of habitable conditions can range from thousands to hundreds of mil- 
lions of years, although planets left with exactly the right amount of hydrogen 
for permanent habitable conditions only occur in rare cases. 

Several effects that we have not considered could alter the simple analysis 
presented here. For one, we have assumed a short, discrete phase of planetary 
accretion, but meteorite impacts can continue long after the main stages of 
formation are over. Events such as the Late Heavy Bombardment in the Solar 
System could erode primordial hydrogen atmospheres and evaporate oceans, 



potentially destroying any emerging life. However, as noted by |Abramov and 



Mojzsis| [2009 , even on an Earth-mass planet relatively large impactors (r > 



50 km) are incapable of sterilizing the entire surface. 

Perhaps more restrictive is the fact that we have only described the evolu- 
tion of pure H2 atmospheres. For more realistic hydrogen-dominated atmo- 
spheres containing a mixture of gases, other sources of opacity will be present 
that could increase greenhouse warming. In addition, photochemistry in the 
planet's upper atmosphere can lead to the formation of aerosols and organic 
hazes, which can increase the planetary albedo or cause 'anti-greenhouse' ef- 



fects [e.g., McKay et al. 1991 . These processes could quantitatively alter 
the results described here, and should be investigated in future using more 
detailed models. Nonetheless, it is difficult to see how they could invalidate 
the basic mechanism entirely. The effectiveness of hydrogen as a greenhouse 
gas means that planets with primordial H 2 envelopes should have high initial 
surface temperatures even if they form far from the host star. In cases where 
the hydrogen envelope is later removed and the planets receive significantly 
less starlight than Earth, transient periods of habitable surface temperatures 
of some duration should therefore always result. 

The reducing chemistry of hydrogen-rich atmospheres provides fertile con- 
ditions for the formation of some of the basic building blocks of life, such as 
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amino acids and ribonucleotides Miller, 1953, Ferris et al. , 1978, Powner et al. 



2009 . Provided C, N and O-bearing species are also present, an exoplanet with 



a slowly evaporating hydrogen atmosphere would likely therefore experience a 
steady rain of organic compounds towards the surface due to photolysis in the 
upper atmosphere. If the atmosphere were dense enough to cause high tem- 
peratures at the surface, these compounds would be destroyed by thermolysis 
relatively rapidly. However, as the greenhouse effect of hydrogen diminished, 
temperatures would decrease, allowing them to begin to accumulating on the 
surface. Moderately high temperatures (60 - 100 °C) dramatically accelerate 
key reactions in primordial biochemistry jStockbridge et aL 2010 , so the ini- 
tial stages of a planet's transiently habitable period might provide particularly 
favourable conditions for the development of complex reaction pathways on 
the surface. 



In an environment dominated by hydrogen-bearing species, the efficiency of 
redox reactions available for metabolism would clearly be very different from 
that of either present-day Earth or the weakly reducing N 2 -C02-H 2 atmo- 
spheres generally believed to be representative of the Archean [e.g. Kasting 



1993 . Nonetheless, emergent organisms could potentially extract energy from 



a variety of mechanisms. Volcanic emission of gases such as CO2 or sulphate 
aerosols from an oxidising mantle Elkins-Tanton and Seager, 2008 would al- 



low for some of the redox reactions already known to be used by methanogens 
and sulphur-reducing bacteria on Earth. Photochemistry itself provides fur- 
ther possibilities: by definition it is driven by the free energy from photons and 
as such results in compounds that are out of equilibrium with their environ- 
ment Pascal and Boiteau 2011 . To give one example, dissociation of water 



vapour and the subsequent loss of hydrogen to space results in highly reactive 
radicals such as OH*, which can cause the formation of oxidized trace species 
(e.g. H 2 2 ) that are transported to the surface via rainout. 



After the loss of the hydrogen envelope in these scenarios, heavier gases such 
as methane or carbon dioxide could remain in the atmosphere, depending 
on their volatility and escape rates. If these gases provided insufficient green- 
house warming to keep temperatures above °C, surface liquid water would no 
longer be possible globally, although processes such as vulcanism and geother- 
mal heating would still allow locally habitable regions and hence some further 
opportunity for life to develop. Investigating the end state of atmospheric 
evolution self-consistently in these cases, while challenging, would be an inter- 
esting subject for future research. 



If life readily emerges under these conditions but cannot adapt to the changing 
climate by altering either itself or its environment, transient hydrogen green- 
house warming may be of little relevance to the search for life around other 
stars. In this case, constraints based on abiotic climate stabilization mecha- 
nisms such as the carbonate-silicate cycle Walker et al. , 1981 will probably 
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determine the frequency of life-supporting planets. The carbonate-silicate cy- 
cle requires plate tectonics for CO2 recycling, and the extent to which this 
phenomenon depends on a planet's mass, internal heat sources and man- 



tle composition is still under investigation Valencia et al. 2007, Korenaga 



2010 . If only a sub-set of rocky planets in the CO2 habitable zone have ac- 



tive carbonate-silicate cycles, the prevalence of life-supporting planets in the 
stellar neighborhood could therefore be small. 



If, however, the adaptability of life to diverse habitats through evolution and 
environmental modification is high, the frequency of conditions allowing its 
formation may be the critical factor. Life on Earth is adapted to a wide range 
of physical and chemical environments, and has colonized all regions of the 
planet where even minimal liquid water is present, including the Antarctic sub- 



surface Cary et al. 


2010 


. It has been shown elsewhere Laughlin and Adams 


2000 


Abbot and Switzer 


2011 that even planets receiving essentially no stel- 



lar flux are able to maintain ice-covered liquid oceans through geothermal 
heating alone. Hence even in the absence of other warming processes, local- 
ized habitable environments would remain in the time after transient hydrogen 
warming had ceased. 



Since the Archean era, life has also extensively modified the composition of 
the atmosphere. The extent to which this modification has regulated surface 
temperatures to optimize planetary habitability remains a subject of intense 
debate [e.g., Lovelock and Margulis, 1974, Kirchner, 2003 , and resolving the 



argument has proved difficult, in part due to our inability to study the Earth's 
atmosphere in different epochs directly. However, it has been shown here that 
appropriate conditions for the origin of life may occur in a fairly wide range 
of cases. Hence if living ecosystems are capable of regulating global climate, 
even planets outside the standard COvwarming habitable zone might retain 
biospheres for longer periods. The boundaries of the CO2 habitable zone are 
fairly well-defined as a function of star type, and global biospheres may be 
detectable from e.g., the presence of non-equilibrium gases in a planet's atmo- 
sphere. Given the planned future searches for nearby planets with biosigna- 
tures, therefore, in theory it could eventually be possible to assess the relative 
importance of abiotic vs. biotic climate stabilization processes directly. 
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Fig. 1. Example temperature profiles used in the climate calculations for pure H2 at- 
mospheres (solid lines) compared with those for mixed H2-H2O atmospheres (dashed 
lines) . 
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Fig. 2. Calculated effective OLR (F out /aT%) function of total H2 pressure for 
pure H2 (solid line) and mixed H2-H2O (dashed line) atmospheres. As can be seen, 
significant differences between the two cases only occur below around pn 2 = 0-5 bar. 
The empirical fit to the pure H2 curve given by ([!]) is shown as a series of crosses 
between 0.5 and 10 bar. 
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Fig. 3. Total amount of hydrogen lost due to XUV energy-limited / Jeans thermal 
escape as a function of planet mass and orbital distance a) in a G-class system after 
100 MYr and b) in an M-class system after 1 GYr. For reference, the dotted red 
line shows the critical cooling limit of Koskinen et al. 
planets under saturated XUV conditions. 
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Fig. 4. Duration of habitable conditions Ai^aft as a function of planet mass and or- 
bital distance for various starting atmospheric hydrogen inventories, a) and b) show 
results for G- and M-class systems, respectively. Results are not plotted for situ- 
ations where insufficient hydrogen is removed to allow surface temperatures lower 
than T max = 333 K. The dashed blue line shows the outer edge of the COvwarming 
habitable zone for the given stellar luminosity, while in a) the thin solid lines show 
the orbits of Earth and Mars. The black region shows cases for which the atmo- 
sphere is eroded just the right amount to allow permanent surface temperatures in 
the 273 to 333 K range. 
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Fig. 5. Total amount of hydrogen lost due to atmospheric Jeans escape only after 
10 GYr as a function of planet mass, for planets with constant exospheric tempera- 
tures. For these relatively low values of T exo , efficient hydrogen dissociation was no 
longer assumed, and the Jeans flux ^ was calculated using the collision cross-sec- 
tion for molecular hydrogen <th 2 = 7t( 1.2 x 1(T 10 ) 2 m 2 , to give a simple lower limit 
on the total escape. 
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